Energy level statistics in weakly disordered systems: 
from quantum to diffusive regime 
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We calculate two-point energy level correlation function in weakly disorderd metallic grain with 
taking account of localization corrections to the universal random matrix result. Using supersym- 
metric nonlinear a model and exactly integrating out spatially homogeneous modes, we derive the 
expression valid for arbitrary energy differences from quantum to diffusive regime for the system 
with broken time reversal symmetry. Our result coincides with the one obtained by Andreev and 
Altshuler [Phys. Rev. Lett. 72, 902 (1995)] where homogeneous modes are perturbatively treated. 
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£ I I. INTRODUCTION 

Much research interest has been attracted by the universal properties of quantum systems with randomness; the 
randomness arises from various sources such as complexity of the system for the case of, e.g., complex nuclei, stochas- 
ticity for disorderd conductors and instability of classical trajectories for chaotic systems [1]. If wave functions have 
enough time to diffuse throughout the system, the spatial strucure of the system becomes immaterial. For such cases, 
I ' the system is called to belong to the ergodic regime where the statistical properties have been extensively studied. For 
example, correlation functions of the energy levels [2] or the scattering matrices [3,4] are shown to have the universal 
forms where specific features of individual systems are represented by a few parameters [1]. These universal forms 
can be calculated with use of standard random matrix models belonging to pertinent Gaussian ensembles. 

If the relevant time scale is comparable with the diffusion time, we need to turn to spatial dependent parts of 
the correlation functions which are not included in the universal random matrix results. These spatial dependent 
CNI . parts for weakly disordered systems are the subject of this paper. The weakly disorderd systems are characterized by 
the condition that j > 1 where g is the dimensionless conductance. The spatial dependent parts of the correlation 
functions are expected to appear as higher order corrections with respect to 1/g to the universal random matrix 
results. The two-point energy level correlation function R(s) (see Eq. (1.1)) with corrections up to 1/g 2 has been 
calculated by several authors: In their seminal work, using diagrammatic perturbation theory, Altshuler and Shklovskii 
[5] obtained the smooth part of R(s) valid for large energy differences s 3> 1. Deviations from the universal statistics 
CNI \ are determined by the diffusive classical dynamics of the system. However, diagrammatic perturbation does not give 
■ an oscillatory part of R(s) characteristic of the universal random matrix result. 

The 1/g 2 corrections including the oscillatory part were first obtained by Kravtsov and Mirlin (KM) [6] using the 
supersymmetric a model [7], the effective theory of a particle moving through a disorderd conductor. Their result is, 
however, valid only for small energy differences s « j [8]. Later, also using the supersymmetric a model, Andreev 
and Altshuler (AA) [9] claimed that the oscillatory part can be obtained from perturbation around a novel saddle 
point in additon to the ordinary one. Their result relies on perturbation and hence is restricted to the region of large 
energy differences s ^> 1. The AA's result was subsequently reproduced with use of the replica [10] and Keldysh 
[11] a model. In quantum chaotic systems, the analysis based on the semiclassical periodic orbit theory [12] and the 
ballistic a model [13] led to similar results to the AA's. 

In all these results, deviations from the universal form is expressed in terms of eigenvalues of a certain operator or its 
spectral determinant. This operator is a diffusion operator for the case of disorderd conductors or a Perron- Frobcnius 
operator for quantum chaotic systems [14]. We can surmise from these results that universal forms also exist in the 
deviations from universal random matrix resuls where parameters representing system specific features are now the 
spectra of this operator. 

At present the available results are yet insufficient because KM's result is valid for s <C g and AA's one is valid for 
s ^> 1. The method to calculate an expression valid for entire range of energy differences is desirable [15]. This paper 
presents one such method. Extending the KM's procedure, we calculate 1/g 2 corrections to the two-point energy level 
correlation function 

R(s) = A 2 (p(E)p(E + u;)) (1.1) 
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valid for arbitrary u>, where s — u)/A, A is the mean level spacing, p{E) is the density of states at energy E, and (. . .) 
denotes averaging over realization of the random potential. We consider disorderd systems with broken time reversal 
symmetry because the pertinent random matrix ensemble (unitary ensemble) is the simplest among the three classical 
ensembles (unitary, orthogonal, and symplectic). We follow the KM's procedure but do not use u expansion which 
limits the range of their result's validity to small u. As a result, we verify that the AA's expression gives the correct 
j > 1 asymptotic at arbitrary to for unitary ensemble. For the unitary case, invoking a mathematical theorem, it 
is implicitly shown [16] that AA's result is exact in spite of their use of perturbation. Here, we want to show it by 
explicit calculation which we expect to be applicable in other ensembles with some modifications. 



II. MODEL AND METHOD 



The supersymmetric a model is a field theory whose variables Q are supermatrices [7]. Statistical properties of a 
particle moving through random potentials can be derived from a generating functional of Q. The expression for R(s) 
reads 



R(s) = -^_Re J dQe~ s ^ J drstikAQ(r) 



(2-1) 



S(Q) = ^ I dfsbr 



^(Vg(r)) 2 +2i S +AQ(f) 



(2.2) 



Here str denotes the supertrace, Q(f) = T(r) -1 AT(r) is a 4 x 4 supermatrix with T(r) belonging to the coset space 
U(l, 1|2), A = diag(l, 1, —1, —1), k = diag(l, —1, 1, —1), V is the system volume, D is the classical diffusion constant, 
s + = s + i0 + , and + is positive infinitesimal. We use the notaions of Ref. [7] everywhere. 

The ordinary saddle point is Q = A and AA's novel one is Q = — fcA. The expansion around these saddle points and 
calculation of the functional integral by Gaussian approximation gives AA's result. The quadratic term of the action 
S(Q) defines the diffusion propagator in terms of which a pertubation series is constructed by Wick's contraction. 
When s goes to zero, the diffusion propagator diverges for spatially homogeneous modes (which we hereafter call the 
zero-momentum modes, or zero modes for short, because their momentum is zero). Thus their approach is restricted 
to the region s ^> 1. 

KM have overcome the dfficulty of zero mode divergence by treating the zero modes and the other spatially 
inhomogeneous modes with non-zero momenta (heareafter called non-zero modes) separately. The matrix Q(r) is 
decomposed in the follwing way [6]: 

O(^=r - 1 O(r)T , (2.3) 

where T is a spatially homogeneous matrix from the coset space and Q(r) describes all non-zero modes (see Eqs. (2.4), 
(2.5)). After perturbative calculation of the contributions by non-zero modes, there remains a defnite integral over a 
few zero mode variables, which can be evaluated without use of perturbation. This amalgamation of both perturbative 
and non-perturbative calculation, in principle, could give the weak localization corrections valid for entire range of s. 
In the course of calculation, however, assuming small u>, KM expand the exponential of the energy term (the second 
term of the action (2.2)). Thus their result is restricted to s <C g. 

The difficulty arising when the energy term is kept on the shoulder of exponent is that the propagators of non- 
zero modes become dependent on zero mode variables. Then after eliminating non-zero modes by perturbation, the 
resulting zero mode integral contains infinite sums and products of non-zero mode propagators (see Eq.(5.5)) as an 
integrand. Hence the zero mode integral, though still being defnite over a few variables, does not seem feasible. In this 
paper, we show this is not the case: by replacing the order of integration, it turns out that we can first integrate out 
the zero mode exactly. Then integration over the non-zero modes by perturbation can be done without any trouble. 

We resume a summary of notations: When A<i? c (where E c = D/L 2 and L is the system size), the matrix Q(r) 
fluctuates only weakly around A. Then the fluctuations can be treated by perturbation with the expansion parameter 
1/g (where g = E c /A). A convinient parametrization is 

Q(r) = (1 - W(r)/2)A(1 - W{r)/2)~\ (2.4) 

where 
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W2i(r) — kWi2(r)^ . W^if) is expanded as W\2{r) — 4>q(?)W\2^ where <j>^ is an eigenfunction of the diffusion 

operator — Z?V 2 with an eigenvalue D\q\ 2 . {4>q{r)} constitute a complete orthonormal set and (j>o{r) = \j\f~V . Thus 
/ dfW(r) oc J drW(r)(f>o(r) = 0. The Jacobian of the transformation Eqs. (2.3) and (2.4) from the variable Q to 
{T ,W} is 

J(W) = 1- — drd? (sir W 12 {r)W2i{r' )f + 0{W b ). (2.6) 

The derivation is given in Appendix A. 

The spatially homogeneous supermatrix Tq is parametrized in a quasi-diagonalized form 

T = U-'foU (2.7) 

with 'eigenvalue' matrix T given by 

~ / cos | isinfe 1( M , , 

ro= ^sta«e-i* cos* J' (2 ' 8) 

^( 9 ,fl). #-K' Q _l 



o i0 B y ' ^ v / ' 

where < 9b < oo, < #f < tt, < tps < 27r, < ipp < 27r. The 'diagonalizing' matrix [/ is given by 

"l = cxp ^ ^ ^ ^ , v 2 = cxpi ^ ^ ^ 2 ^ , (2.11) 

where £i and £2 are anticommuting variables. For this paramctrization, the measure effo is given by 

dT = — smh °b sin 6 F ddBdeFdlf>BdipFd ^* d ^*_ (2.12) 
(coshes — cos9p) 



III. CHANGING THE VARIABLES 

To simplify the zero mode integration, we change the variables. Substituting the decomposition (2.3) into Eq. (2.1), 
we obtain 

str (VQ(r)) 2 = str (vQ(r))' , (3.1) 
J drstiAQ(r) =strAT - 1 YT a , (3.2) 

J dr str fcAQ(r) = str kATo^YTo, (3.3) 

where Y = J dfQ(r). We make the supermatrix Y a block diagonal form as follows (for detail see Appendix B): 

e x Ye- x = Q, (3.4) 

where A is a block off-diagonal supermatrix and 
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Q = A (v + 1 y dr TU 2 (r) + i ^ drW 4 (r) + 0(TU 6 )^) . 



(3.5) 



Up to 1/g 2 order calculated here, e x T can be replaced with T (see Appendix B). Thus in Eqs. (3.2) and (3.3), Y 
can be replaced with Q. 

To eliminate the anticommuting variables of the zero modes in the action S(Q), we change the variables as 
UWU- 1 -> W. Eventually, we get 



J df str AQ(f) = str Af -1 QTb, 



(3.6) 



/ 



df str kAQ(r) = str kAU^T^QT^U. 



(3.7) 



IV. INTEGRATION OVER THE ZERO-MOMENTUM MODE VARIABLES 



Expanding the pre-exponential term of Eq. (2.1) in the anticommuting variables of the zero modes, we obtain 



(strMC/- 1 f - 1 Qf f/) 2 - (strfcAf-^fo) 2 - 2£6£& (str hf^QToj 



+ 



(4.1) 



where we have used that strQ = and the dots indicate terms which vanish after integration of the zero mode 
anticommuting variables. We write the contribution by the first (second) term of Eq. (4.1) to R(s) as R\(s) (^(s)). 
With use of Parisi-Sourlas-Efetov-Wegner theorem [17], integration of zero mode variables in R\(s) becomes the value 
of its integrand at 9 B = 9 F = 0: 



Ri(s) = ^72 Re J dWJ{W)e~ s ^ (strfcAg) 2 + 0(l/ ff 3 ), 



where 



S(W) = -str 



^Jdr (VQ(0) 2 + 2Ls+ AQ 



(4.2) 



(4.3) 



Integrating the zero mode anticommuting variables, we obtain 

R 2 ( S ) = -_L Re f dW d9 B de F J (W) sinhfl B8 infl F e _ s{foW) U A f-iQ f \ 2 + 0{l/g% (<u) 
ov J (cosh 6b - cos Of) v ' 



S(T ,W) = — str 



- J df (vQ(f)) 2 + 2is+Af - 1 QT 



Here. 



where P(B) 
strP(B)AQ -- 



where 



str Af^QTo = cosh 6 B str P(B)AQ + cos 6 F str P(F)AQ, 



(4.5) 



(4.6) 



(1 — k)/2 and -P(i^) = (1 + fc)/2 and 1 is an unit matrix. Using coshes = Xb, cosdp = Xp, 
-2V (1 + A), and strP(F)AQ = 2V(1 + B), we obtain 



i? 2 (s) = -Re / dW J(W)exp 



nD 
4UA 



y drstr (VQ(r))' 



7( S ,TU)+0(l/ 5 3 ), 



/OO /■ 1 

dA B y d\ F f(\ B ,\ F )g(\B,\F), 



(4.7) 



(4.8) 
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/(A B ,A F ) 



J iirs+(\ B -\ F ) 



(Xb-XfY 



g(X B , X F ) = e^ s+ ^ B a—\f b) [Xg (1 + A) _ Xp (1 + B y2 
We introduce an indefinite integral of f(X Bl X F ): q\%\ b F(Xb, X f ) = f(X B ,X F ), 

/Xf poo . i 
db / dae 1 ^^- — 
-oo Jx B (a-b) 2 

/OO 
dx e i,TS+ ( 



(4.9) 



(4.10) 



1 1 



(4.11) 



where we have used the transformation (a — b) 2 = / °° te ^ a b ^dt and i7rs+ — t = i7rs+x. Then partial integration 
gives 



I(s,W)= [F(Xb,Xf)9(Xb,Xf)} X x b b ZT 



A F = 1 



A F = -1 



d 

dX B F(X B , X F )——g(X B ,X F ) 
o\ B 



A B =oc 



X F =1 

Af=-1 

d 2 



dX F F(X B ,X F )—g(X B ,X F ) + dX B dX F F(X B ,X F )———g(X B ,X F ). (4.12) 

-1 OA F \\ B =l Jl J-l OA B 0A F 



After A-integration, we obtain 



I(s,W) = - 



Atts+(Xb A-X f B) 



A b = 1 



A F = 1 



\T!S+(A-B) 



dxe ins+ ^ B - x ^ x l(X B ,XF,x) 

/OO p OO 

dxl{l,l,x) ^ e ^ s+ ( A + B ) J dxe 2l7TS+x l(l,-l,x), 



X F = -1 



(4.13) 



where 



I(X B , \ F ,x)=(--\) A- Y B- X [X B (1 + A) - X F (1 + B)f 



x x" 



2 n.i u -^ 



i7r,s+ V a; 2 x 3 



A(l + B) + S(l + A) [A B (1 + A) - X F (1 + B)\ 



+ 



1 L ) i- 3 i?- 3 

(i7r,s+) 2 \x 3 x 4 



A 2 (1 + B) + AB (1 + A) (! + £?) + B 2 (1 + A)' 



(4.14) 



Here, A = 1 + ^ and _B = 1 + -f • From Eq. (4.13) we can read that the propagator do not include zero mode variables. 
Therefore the difficulty of KM's approach mentioned in Sec. II is resolved. We write the contribution by the first 
(second) term of Eq. (4.13) to R 2 (s) as i?^ 14) (s) 



V. INTEGRATION OVER THE NON-ZERO-MOMENTUM MODE VARIABLES 

In order to calculate the remaining integral over the non-zero mode variables W with use of the Wick theorem, we 
need to know the contraction rules for the action 



So = — I df str 



D 



(W(r)) 2 + is+ (X B P(B) + X F P(F)) W(ff 



The result is summarized as 

(stTXW(r)yW(r>)) = -i ]T II(r, r%, [str P(g)X str P(g')y - str AP{g)X strAP{g')y] , 



(5.1) 



(5.2) 
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(strAW(r) strJW^)) = -J^n(f,rV [str P(g)XP(g')y - str AP(g)XAP(g')y] , (5.3) 



where X and y are arbitrary 4x4 supermatrices, g and g' are £> or F, 

(...) =P- 1 ( S ,A B ,A F ) J dW(...)e- 3 ", 

<nr \ \ \ tt n(g; A B , A B )n(<7; A F , A f ) 
7 • A ' i - A " H n (? ;A B ,A F) 2 



(5.4) 
(5.5) 



n (9 ;A g)V ) = ^(^-i S ^ 



(5.6) 



n(r, = ^^-(f)^(rOn(g; A s , A s ,)- 

For i?i(s), using the contraction rules with A# = 1 and Af = 1 we obtain 

2\ 



= 1 + ^ Rc E n ^ 1 ' 1 ) 2 + °( 1 /s 3 ) 



o(i/ ff 3 ) 



(A B ,A F ) = (1,1) 



9#0 



l-^^lnP( S , 1,-1) + 0(1/^). 



(5.7) 



(5.8) 



In order to get i?2(s), we need to estimate the order of the each term of Eq. (4.14) with use of the fact that (i) 
W 2 <~ 0(1/ g) for all s, (ii) the property of the exponential integral 



E, 



-dx 



where Re a > and n e N, and (hi) KM's result [6]: for s < 1, 



where C n (s) is 0(1). 



n=0 y 



A. (A s , Af) = (1, 1) point 



(5.9) 



(5.10) 



Expanding 1(1, 1, x) in W, there appears a term proportional to 1/x, (A — B) 2 /x. The x-integration of this term, 
at first glance, seems to diverge. Actually it causes no problem becaue it vanishes after W-integration. In order to 
obtain the corrections up to 1/g 2 , it is sufficient to expand the remaing terms up to the 4th order of W for entire 
range of s [18]; 

1(1,1, X ) = -±(A-B) 2 

2 (4-4) ( A - B )+( 1 2-^ + ^) ( A2 - B2 ) 
\x x° J \ r x x I 



(ilTS+) 2 

+0(l/g 3 ). 



1 3 2 

x 3 x 4 - X 1 



1 9 18 10 

x 3 x 4 - x 5 x' 



(5.11) 



G 



After x-integration, we obtain 

/ 1(1, l,x)dx = -(A Bf —-(A - B) + —— + 0(1/ g 3 ). 
After IF-integration, we obtain 

4 M) (*) 



i 



2tT 2 S 2 ' 



(5.12) 



(5.13) 



B. (\ b ,\f) = (1,-1) point 

For s > 1, E n (— i27rs+) in Eq. (4.13) is 0(1). Accordingly, by the same reason as the case of (Xb,^f) = (1,1) 
point, it is sufficient to expand 1(1, — l,x) up to the 4th order of W for entire range of s again. 



l(l,-l,x)=4 
+4 
+ 



1 1 \ 4/1 



1 



6 / 1 



1 





2 






G- 









2 5 3 

~2 3 4 



1 5 7 

+ 



3 \ 4/1 



(i /' I 3.2 
6 / 1 



1 1 



4 / 1 



12 6 

x 

3 2 \ 2/1 9 18 10 



9 18 10 

+ 



[A -Bf 



+0(l/g 3 ) 

where s = —i2ns + . Using the relation E n (s) + jE n+ \(s) = e~ s /s, we obtain 

1 3 



dxe-' sx l(l,-l,x) = 4^ + 1 



2 I dxe- sx 



(A + Bf 



(A-B) 2 + 0(l/g 3 ), 



(5.14) 



(5.15) 



POO IZTTS 

\,e+(A+B) / dxe i2ns+ XI ^_^ x)= i,s+(A^ + B^)e 

Jl (17TS+ 

+ l7TS+(A^+B^>) 

+0(l/g 3 ), 



127TS + 



1TTS 



(5.16) 

where T(s) = e [27TS+ + 2 dxe l2lTS+x - Jr). The contribution by the second term of Eq. (5.16) is 0{l/g 3 ), 
because for s > g it is obvious and for s < g it is found after lU-integration. 

By the same way, the contributions by 0(W i ) terms emerging from ^VQ(r)j and the jacobian J(W) turn out to 

be 0(1/ g 3 ). Finally, we obtain 



R i 2 1 '~ 1 \s) = --Rc I dWJ(W)exp 



irD 
4UA 



J rfrstr (V0(r)y 



(-+X2 + 



cos 2-ks 
2tt 2 s 2 



V(s,l,-l) + 0(l/g 3 ). 



(5.17) 



VI. RESULT AND SUMMARY 



Equations (5.8), (5.13), and (5.17) give 
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= i - + ^S(-) + o(i/*»), (6.i) 

where X>(s) = T>(s, 1, — l)/s 2 is the spectral determinant of the classical diffusion operator, in agreement with 
Refs. [9-11]. Here the contribution of the (As, Xf) = (1, 1) 

-i^ ln ^)=8^ Re E n ^ 1 ' 1 ) 2 ( 6 - 2 ) 

<? 

is the ordinary perturbative result [5] . 

In summary, using KM's separate treatment of the zero and non-zero modes, performing the integration over the 
zero modes exactly and subsequent integration over the non-zero modes perturbatively, we find the expression for the 
two-point energy level correlation function valid for arbitrary energy differences up to 1/g 2 order for unitary ensemble. 
This expression is same as the one [9-11] obtained in a saddle-point approximation which in general valid only for 
u> 3> A [19]. By explicit calculation, we verify that the exactness of the saddle-point answer for unitary ensemble, 
which is guaranteed in the ergodic regime due to the Duistermaat-Heckman theorem [16], holds even in the diffusive 
regime. Since this is a specific feature of the unitary ensemble, it is very interesting to investigate the expressions for 
the other (orthogonal, symplectic, etc.) ones. 

The reason why AA's novel saddle point appears is as follows: The zero-mode integral is carried out exactly. Then 
the remaining expression for non-zero modes is evaluated at terminal points of the integral region, 1 < < oo, — 1 < 
Xf < 1- These terminal points precisely correspond to the ordinary (Xb = Xf = 1) and AA's novel saddle point 
(Xb = 1, Xf = —1). This correspondence between the terminal points of zero-mode variables and saddle points seems 
to hold for the orthogonal ensemble. Hence, for the orthogonal case, a similar way of calculation may work: although 
exact integration of zero mode variables is no longer probable, the iteration of partial integration enables one to 
evaluate the integral on the terminal points up to a necessary order of 1/g, if the additional terms generated by these 
iterations are only higher order term of 1/g. For the symplectic case, the correspondence seems more subtle because 
the novel saddle points are not isolated but consist of a saddle point maifold. Further investigation will be necessary 
for the symplectic case. 
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APPENDIX A: CALCULATION OF THE JACOBIAN 



We give the calculation of the Jacobian J of the transformation Q — ► {W, T } because, to the best of our knowledge, 
there is no explicit expression of it in literature. It is calculated by the following Gaussian integral : 



J d(SW)d(8T ) exp - J drstr (SQ(f) f 

Here SQ(r) is variation of Eq. (2.3) : 

5Q(r) = To 1 (§Q(r) + [o(r),^ 

where ST^ = STqT^ 1 . Then, 



J- 1 . 



(Al) 



(A2) 



J drsti {5Q{f)f = j drstr {sQ{r)^ + 2 j drstr5Q(r) Q(r),6T^ + J drstr Q(r),ST^ 2 . (A3) 



The integrand of Eq. (Al) is expanded up to 4th order in W. 
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exp 
= exp 



The 5W dependent part of the integrand is 

-J dfstr (SQ(r)Y -2 J dfstr 5Q(f) Q(r),ST^ 
J drstr (SW) 2 

x|l + iy dfstrW 2 {5W) 2 + 2 drstr X 1 (SW)ST^ > 

+2 (^j dfstr X 1 (SW)ST^j (^j dfstr X 2 (SW)ST^j 

+ 1 ( drstr W 4 (SW) 2 + \- [ drstr (W 2 SW) 2 + I 
8 J 16 J 8 

dfstr X 1 (SW)ST^j (^j dfstr X 3 (SW)ST^j 



dfstr W 2 (SW) 2 



~ (J dfstr X 2 (SW)ST^j df str X 1 (SW)S% 



dfstr W 2 {SWf 



(^J dfstr X 1 (SW)ST^j +0(W 5 



where 



X n (SW) = Y^{-^) k W k 5WW n - k . 



k=0 



The SW independent part is 



exp 



drstr 



= exp [2^str (6T^ - AST^A)ST^\ 

x |l - 2 J dfstr A<5TqA (W 2 ST^ - WSTqW) 

- J dfstr AST^A (W 3 8T^ - 2W 2 ST' Q W) 

_I J dfstr ASTqA (W i 5T[ ) - 2W 3 ST^W + W 2 ST^W 2 ) 



+2 



J dfstr ASTqA (W 2 5T[ ) - WST^W) 



0(W 5 ) 



(A4) 



(A5) 



(A6) 



The integral of the SW dependent part over SW can be calculated by using the Wick theorem and the contraction 
rules: 



1/1 



{str X5W{f)ySW{r')) sw = - I — - S(f- r') ) (str X str y - str AX strA^) : 



(A7) 



1 / 1 



(str XSW(f) str ySW(r'))sw = ^ I y - S(f - r') ) (str Xy - str AATA?) . 



where X and y are arbitrary supermatrices and 

(. . .) 6W = J d(SW) (. . .) exp J dfstr (SW) 2 



(A8) 



(A9) 



+- 



+ 



The result is 

J d(SW) exp - J drsti (<SQ(f)) 2 -2 J drsti SQ(r) Q{r),5T^ 
= 1 + J drstr {ST^ + ASTqA) (W 2 6T^ - WST^W) + j dr str AST^A (W 3 ST^ - 2W 2 8T^W) 
_I j drsti (W 2 ST[ ) ) 2 + ^ j dr str AST^A (4W A ST^ -8W 3 5T[ ) W + 5W 2 6T[ ) W 2 ) 

j dr str {5T^ + AST^A) (W 2 8T^ - WST^W) 

7 J dfd? str (ST^ - ASTqA) (2W {rfW {r'f 5T^ - 2W{r) 2 W{P)8T^W{?) 
-2W(r)W(r') 2 ST^W(P) + 2W(r) 2 5T^W(P) 2 + W{r)W{P)5T^W{P)W(r)^ 

-^V J d r(^ W2 ) 2 + 64^(7 drstiW 2 ^ + 0{W% 
Multiplying this by Eq. (A6), we find 

J d(5W) exp - J drsti (SQ(r)) 2 
= exp [21/str (5T£ - AST^ST^} 

x |l + J dr str (STo - AST^A) (W 2 ST^ - WST^W) 

+ i j dr str (5T^ - AST^A) (\V 2 ST^ - WST^W) + dr str W 2 ^j 

+ W J d?dr: ' sil ( ST ° ~ A5T o A ) (2W(r) 2 W{P) 2 5T^ - 2W{r) 2 W(P)ST^W{P) 
-2W(r)W(P) 2 ST^W(P) + 2W(r) 2 ST^W(P) 2 + W(f)W(P)6T^W(P)W(rfj 
_I J dr str (ST^ - ASTqA) W 2 ST^W 2 - J dr (str W 2 f + 0(W 5 )^ . 

The remaining ST integral can be also calculated with use of the Wick theorem and the contraction rules: 

(str X(ST^ 12 y(ST^ 21 ) S To = ~sbr X str?, 



(A10) 



(All) 



(A12) 



(str*(J22)i2 stry(<5^) 2 i) 4 r = -g^str*?, 
where X and y are arbitrary supcrmatrices and 

(. • Vo = / d(ST ) (. . .) exp [2Vstr (ST^ - AST^A)ST^] . 



(A13) 



(A14) 



Since the zero mode variables enter into Eq. (All) only through (5Tq) 12 and (STq)i 2 , the Jacobian does not contain 
the zero mode variables. We finally obtain the Jacobian (2.6) from the expression 



/ 



d(SW)d(ST Q ) exp 



rffstr (8Q(r)Y 



= 1 



16V 2 



J drdP (str W 12 (r)W 2 i(P)y + 0(W 5 ). (A15) 
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APPENDIX B: BLOCK DIAGONALIZATION 



We can diagonalizc Y = J drQ(r) in the retarded-advanced space as follows. For Eq. (3.4), 

Y = ^Y {n \ (Bl) 



n=0 



X = Y J * (2n+1 \ {x (2n+1 \A} = 0, (B2) 



n=0 



Q = Y,Q i2n \ [0 (2n) ,A]=0, (B3) 



n=0 



where Y {n \ X {n \ Q (n) mean 0(W n ) term for Y, X, Q and {, } stands for an anticommutator. Then Q and X are 
perturbatively determined and given by Eq. (3.5) and 

X = ^ J drW 3 {r) + 0{W 5 ). (B4) 

We can parametlize e x T in terms of T and R: 

e x T = TR, (B5) 

where T and R obey the properties 

T -1 = ATA, [R,A}=0, (B6) 

and are given by 

T = (e x T*e x ) 1/2 = To + 0(W 3 ) + 0(W 5 ), (B7) 
R = T~ 1 e x (e- x T 2 e- x ) 1/2 = 1 + 0(W 3 ) + 0(W 5 ). (B8) 

From Eq. (B7), 

dT = (l + 0(W 3 ) + 0(W 5 )) dT . (B9) 
Substituting the transformation (3.4) and (B5) in Eqs. (3.2) and (3.3), we find 

J drstiAQ(r) = strAT^QT, (BIO) 

(^J drstrkAQ(r^j = (str kAT^QT} 2 + 0(W 3 ) + 0(W 5 ). (Bll) 

In Eqs. (B9) and (Bll), 0(W 3 ) terms contribute to 0(W 6 ) in R(s) and are negligible. Therefore we replace T with 
T . 



[1] T. Guhr, A. Miiller-Groeling, and H. A. Weidenmiiller, Phys. Rep. 299, 189 (1998). 
[2] E. Brezin and S. Hikami, Phys. Rev. E 56, 264 (1997). 

[3] G. Hackenbroich and H. A. Weidenmiiller, Phys. Rev. Lett. 74, 4118 (1995). 



11 



[4] N. Mae and S. Iida, Phys. Rev. E 63, 047102 (2001). 

[5] B. L. Altshuler and B. I. Shklovskii, Zh. Eksp. Teor. Fiz. 91, 220 (1986) [Sov. Phys. JETP 64, 127 (1986)]. 
[6] V. E. Kravtsov and A. D. Mirlin, Pis'ma Zh. Eksp. Teor. Fiz. 60, 645 (1994) [JETP Lett. 60, 656 (1994)]. 
[7] K. B. Efetov, Super symmetry in Disorder and Chaos (Cambridge University Press, Cambridge, 1997). 
[8] For a random matrix model describing crossover between Poisson and CUE statistics, relevant result to the KM's is 
obtained, see T. Guhr and A. Miiller-Groeling, J. Math. Phys. 38, 1870 (1997); K. M. Frahm, T. Guhr, and A. Miiller- 
Groeling, Ann. Phys. 270, 292 (1998). 
[9] A. V. Andreev and B. L. Altshuler, Phys. Rev. Lett. 75, 902 (1995); A. V. Andreev, B. D. Simons, and B. L. Altshuler, 
J. Math. Phys. 37, 4968 (1996). 
[10] A. Kamenev and M. Mezard, Phys. Rev. B 60, 3944 (1999). 
[11] A. Altland and A. Kamenev, Phys. Rev. Lett. 85, 5615 (2000). 
[12] E. B. Bogomolny and J. P. Keating, Phys. Rev. Lett. 77, 1472 (1996). 

[13] A. V. Andreev, O. Agam, B. D. Simons, and B. L. Altshuler, Phys. Rev. Lett. 76, 3947 (1996); Nucl. Phys. B 482, 536 
(1996). 

[14] O. Agam, B. L. Altshuler, and A. V. Andreev, Phys. Rev. Lett. 75, 4389 (1995). 

[15] The supersymmetric a model description is valid under the condition: u<1/t where r is the elastic mean free time (see 

Rel. [7]). The energy difference u> considered in the paper is assumed to satisfy this inequality. 
[16] M. R. Zirnbauer, cond-mat/9903338. 

[17] M. R. Zirnbauer, Nucl. Phys. B 265, 375 (1986); F. Constantinescu and H. F. de Groote, J. Math. Phys. 30, 981 (1989). 
[18] Eq. (4.14) contains 1/s and 1/s 2 factors. One may wonder what order of 1/g should be assigned to these factors because 

they can become arbitrary large when s goes to 0. However KM's result (5.10) shows these factors combined with all other 

s dependent terms give no divergence at s — 0. Hence we need not take these factors into account for order estimation. 
[19] In AA's calculation, technically speaking, they had to introduce an external parameter in order to avoid the divergence of 

their integral over zero mode variables of the Fermion-Fermion block for the novel saddle point. However, such an artificial 

regulator is not necessary in our calculation. 



12 



